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We consider the neutralino proton cross section for detection of Milky Way dark 
matter for a number of supergravity models with gauge unification at the GUT 
scale: models with universal soft breaking (mSUGRA), models with nonuniversal 
soft breaking, and string inspired D-brane models. The parameter space examined 
includes nij^ <1 TeV and tan/3 < 50, and the recent Higgs bound of >114 
GeV is imposed. (For grand unified models, this bound is to be imposed for all 
tan /3.) All coannihilation effects are included as well as the recent NLO corrections 
to b — » S7 for large tan/3, and coannihilation effects are shown to be sensitive to 
Aq for large tan/3. In all models, current detectors are sampling parts of the 
paramater space i. e. tan /3 ~ 25 for mSUGRA, tan /3 ~ 7 for nonuniversal models, 
and tan /3 ~ 20 for D-brane models. Future detectors should be able to cover 
almost the full parameter space for fi > 0. For fi < 0, cancellations can occur for 

m i/2 ~ 450 GeV, allowing the cross sections to become ~ 10 _1<) pb for limited 
ranges of tan/3. (The positions of these cancellations are seen to be sensitive to 
the value of cr^tf.) In this case, the gluino and squarks lie above 1 TeV, but still 
should be accessible to the LHC if mi/% < 1 TeV. 

1 Introduction 

The existance of dark matter, which makes up about 30% of all the matter 
and energy in the universe, is well documented astronomically. However, what 
it is made of is unknown, and there have been many theoretical suggestions: 
wimps, axions, machos, etc. The Milky way consists of perhaps 90% dark 
matter, and so is a convenient "laboratory" for the study of dark matter, par- 
ticularly by direct detection by terrestial detectors. We consider here the case 
of supersymmetric (SUSY) wimp dark matter and its detection by scatter- 
ing by nuclear targets. In SUSY models with R-parity invariance, the wimp 
is almost always the lightest neutralino, Xii an( A for heavy nuclei, the spin 
independent scattering dominates the cross section. Since then neutron and 
proton cross sections in the nuclei are nearly equal, it is possible to extract the 
X? — P cross section, a^.o_ p , from any data (subject, of course, to astronomical 
uncertainties). Current detectors (DAMA, CDMS, UKDMC) are sensitive to 
cross sections 

CT£o_ p ~ 1 x 10~ 6 pb (1) 
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with perhaps an improvement on this by one or two orders of magnitude in 
the near future. More long range, future detectors (GENIUS, Cryoarray) plan 
on a significant increase in sensitivity, i.e. down to 

<t^o_ p £ (1CT 9 - l(T 10 )pb (2) 

We discuss here how such sensitivities might relate to what is expected from 
supersymmetry models. 

We consider here three SUSY models based on grand unification of the 
gauge coupling constants at the GUT scale of Mg ~2x 10 16 GeV: 

1. Minimal Supergravity GUT Models (mSUGRA)S. Here there are univer- 
sal soft breaking masses occurring at scale Mg. 

2. Non-universal Soft Breaking Modelsi. Here the first two generation of 
squarks and sleptons soft breaking masses are kept universal (to suppress 
flavor changing neutral currents) and the gaugino masses are universal at 
Mg, while nonuniversalities are allowed in the Higgs soft breaking masses 
and the third generation squark and sleptons masses at Mg. 

3. D-brane String Models (based on type IIB Orientifolds)lll. Here the 
SU (2) ^doublet scalar masses are different from the singlet masses at 
Mq, and the gaugino masses are similarly not degenerate. 

The three types of models have varying amount of complexity in the soft 
breaking parameters, and while the first two models arise from natural phe- 
nomenological considerations in supergravity theory, there are also string 
models that can realise such soft breaking patterns. Though physically very 
different, all the models turn out to lead to qualitatively similar results: Cur- 
rent detectors are sensitive to a significant part of the SUSY parameter space, 
and future detectors should be able to cover most of the parameter space ex- 
cept for some special regions where accidental cancellations can occur which 
make o~^o_ p anomalously small. Thus dark matter experiments offer signifi- 
cant tests of supersymmetry over the same time scale (the next 5-10 years) 
that accelerator experiments will. 

While each of the above models contain a number of unknown parameters, 
theories of this type can still make relevant predictions for two reasons: (i) 
they allow for radiative breaking of SU(2) x U(l) at the electroweak scale 
(giving a natural explanation of the Higgs mechanism), and (ii) along with 
calculating cr^o_ p , the theory can calculate the relic density of i.e tt^o = 
p^o/pc where p^o is the relic mass density of Xi an d p c — 3Hq/8ttGn (H is 
the Hubble constant and Gjv is the Newton constant). Both of these greatly 
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restrict the parameter space. In general one has f2^o/i 2 ~ (f f dx(a ann v)) 1 
(where <7 ann is the neutralino annihilation cross section in the early universe, 
v is the relative velocity, Xf = kTf/m^o, Tf is the freeze out temperature, 
(...) means thermal average and h = Ho/100 km s Mpc ). The fact that 
these conditions can naturally be satisfied for reasonable parts of the SUSY 
parameter space represents a significant success of the SUGRA models. 

In the following we will assume Hq = (70 ± 10)km s _1 Mpc _1 and matter 
(m) and baryonic (b) relic densities of Sl m = 0.3 ± 0.1 and f2& = 0.05. Thus 
n^oh 2 = 0.12 ± 0.05. The calculations given below allow for a 2a spread, i.e. 
we take i 

0.02 < fl^oh 2 < 0.25. (3) 

It is clear that accurate determinations of the dark matter relic density 
will greatly strengthen the theoretical predictions, and already, analyses using 
combined data from the CMB, large scale structure, and supernovae data 
suggests that the correct value of the relic .density lies in a relatively narrow 
band in the center of the region of Eq. (3)19. We will here, however, use the 
conservative range given in Eq. (3). 

2 Theoretical Analysis 

In order to get accurate predictions of the maximum and minimum cross 
sections for a given model, it is necessary to include a number of theoretical 
corrections. We list here the main ones: (i) In relating the theory at Mq to 
phenomena at the clcctrowcak scale, the two loop gauge and one loop Yukawa 
renormalization group equations (RGE) are used, iterating to get a consistent 
SUSY spectrum, (ii) QCD RGE corrections are further included below the 
SUSY breaking scale for contributions involving light quarks, (iii) A careful 
analysis of the light Higgs mass rrih is necessary (including two loop and pole 
mass corrections) as the current LEP limits impact sensitively on the relic 
density analysis, (iv) L-R mixing terms arc included in the sfermion (mass) 2 
matrices since they produce important effects for large tan (3 in the third 
generation, (v) One loop corrections arc included to m& and m T which are 
again important for large tan (3. (vi) The experimental bounds on the b — ► 
decay put significant constraints on the SUSY parameter space and theoretical 
calculations here include the leading order (LO) and NLO corrections. We 
have not in the following imposed b — r (or t — b — t) Yukawa unification or 
proton decay constraints as these depend sensitively on unknown post-GUT 
physics. For example, such constraints do not naturally occur in the string 
models where SU(5) (or 50(10)) gauge symmetry is broken by Wilson lines 
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at Mq (even though grand unification of the gauge coupling constants at Mq 
for such string models is still required). 

All of the above corrections are now under theoretical control. In partic- 
ular, the frpjpi sj SUSY NLO corrections for large tan/3 have recently been 
calculated Em. We find that the NLO corrections give significant contribu- 
tions for large tan /3 for fi >0. (We use here Isajet sign conventions for the fi 
parameter.) There have been a number of calculations of o~^o_ p given in the 
literature i@0lll'E10lll, and we find we are in general numerical agreement 
in those regions of parameter space where the authors have taken into account 
the above corrections. 

Accelerator—bounds significantly limit the SUSY parameter space. As 
pointed out in E3, the LEP bounds on the Higgs mass has begun to make a 
significant impact on dark matter analyses-. Since at this time it is unclear 
whether the recently observed LEP events E-3 represent^ Higgs discovery, we 
will use here the current LEP lower bound of 114 GeV t3. There are still some 
remaining errors in the theoretical calculation of the Higgs mass, however, 
as well as uncertainty in the i-quark mass, and so we will conservatively 
assume here for the light Higgs (h) that > 110 GeV for all tan/3. (For 
the MSSM, the Higgs mass constraint is significant only for tan/? ~ 9 (see 
e.g. Igo-Kemenesta) as Ah production with niA = mz can be confused with 
Zh production. However, in GUT models radiative breaking eliminates such 
regions of parameter space and the LEP constraint operates for all tan/3.) 
LEP data also produces a bound on the lightest chargino (xf) of m~± > 102 
V El. For b — + S7 we assume an allowed range of 2a from the CLEO data 

1.8 x 10~ 4 < B(B Y s7 ) < 4.5 x 10~ 4 (4) 

The Tevatron gives a bound of m§ > 270 GeV( for rriq = mg)il. 

Theory allows one to calculate the Xi _ q uar k cross section and we follow 
the analysis of c3 to convert this to x? — P scattering. For this one needs the 
7r — N sigma term, 

o-kN = ^(m u + m d )(p\uu + dd\p), (5) 

Co = 0V at — (m u + md)(p\ss\p) and thCj—.quark mass ratiOj—r = 
m s /(l/2)(m u + m d ). We use here a = 30 MeVE3, and r = 24.4 ± 1.5HL_Rc- 
cent analyses, based on new n — N scattering data gives a^N — 65 MeVt-fEj. 
Older 7r — N data gave ov^r = 45 MeVc3. We will use in most of the analysis 
below the larger number. If the smaller number is used, it would have the 
overall effect in most of the parameter space of reducing x? — P by about a 
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factor of 3. However, in the special situation for /j <0, where there is a can- 
cellation of matrix elements, the choice of ov/y produces a more subtle effect, 
and we will exhibit there results from both values. 

3 mSUGRA model 

We consider first the mSUGRA model where the most complete analysis has 
been done. mSUGRA depends on four parameters and one sign: tuq (universal 
scalar mass at Mq), 7711/2 (universal gaugino mass at Mq), Aq (universal cubic 
soft breaking mass), tan/3 = (H 2 )/(Hi) (where (H( 2 ,i)) gives rise to (up, 
down) quark masses) and Li/\fJ>\ (where /x is the Higgs mixing parameter in 
the superpotential, = fiHiH 2 ). One conventionally restricts the range of 
these parameters by "naturalness" conditions and in the following we assume 
mo < 1 TeV, mi/2 < 1 TeV (corresponding to m§ < 2.5 TeV, m^o < 400 
GeV), |A)/wi/ 2 | < 4, and 2< tan/3 < 50. Large tan/3 is of interest since 
SO(10) models imply tan/3 > 40 and also cr^a_ p increases with tan/3. o"^o_ p 
decreases with rrii/2 for large myi^- 

The maximum value of o~^o_ p arises then for large tan/3 and small m 1 / 2 . 
This can be seen in Fig.l where (a 

v° — T) J max 

is plotted vs. m^o for tan/3=20, 
30, 40 and 50. Fig. 2 shows ^<^h 2 for tan /3 = 30 when the cross section takes 
on its maximum value. Current detectors obeying Eq (1) are then sampling 
the parameter space for large tan/3, small m^o and small Q^oh 2 i.e 

tan/3 ~ 25, m^o ~ 90GeV, Q^oh 2 ~ 0.1 (6) 

Further, as can be seen from Fig. 3, rrih does indeed exceed the current LEP 
bound over this entice region. 

As discussed in til, coannihilation effects in the early universe can signif- 
icantly influence the relic density calculation. To discuss the minimum cross 
section, it is convenient then to consider first m^o ~ 150 GeV (rrii/2 < 350) 
where no coannihilation occurs. The minimum cross section occurs for small 
tan /?. From Fig. 4 one sees 

<7%o_ p ~ 1 x 10~ 9 pb; m^o ~ 140GeV; tan/3 = 6 (7) 

which would be accessible to detectors that are currently being planned (e.g. 
GENIUS). 

For larger m^o, i.e. m 1 / 2 ~ 350 the phenomena of coannihilation can 
occur in the relic density analysis since the light stau, fi, (and also e^, JIr) 
can become degenerate with the Xv The relic density constraint can then be 
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Figure 1. (ffjO_ p )max for mSUGRA obtained by varying Aq and mo over the parameter 

space for tan = 20, 30, 40, and 50 14 . The relic density constraint, Eq. (3) has been 
imposed. 
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Figure 2. fl-oh for mSUGRA when (<r~o _) takes on its maximum value for tan/3 = 30 



satisfied in narrow corridor of mo of width Am-o ~ 25 GeV, the value of mo 
increasing_as mjrt increases and this was examined for low and intermediate 
tan/3 in til. Since mo and mj/2 increase as one progresses up the corridor, 
CvO_„ will generally decrease. 

We consider first the case of fi > E3. Coannihilation effects generally 
begin for mi/2 ~ 400 GeV (m^o ^160 GeV), and it is of interest to see what 
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Figure 3. m/, for mSUGRA as a function of m-o for tan/3 =30, when cr-o_ r , takes on it 
maximum value 14 . 




occurs for large tan/3. For large tan/3, there is only a coannihilation region 
left in the parameter space, and the allowed regions, exhibiting the allowed 
narrow corridors of parameter space are shown in Fig. 5 for tan /3 = 40. In 
this domain the lightest stau (fi) is the lightest slepton due to the large L-R 
mixing in the (mass 2 ) matrix, and so dominates the conanihilation effects. We 
note that the allowed corridors are sensitive to Aq , and large Aq can allow large 
mo as mi/2 increases. The thickness of the allowed corridors also decrease as 
Ao increases. There is also a lower bound on m 1 / 2 for the allowed regions due 
to the b — > S7 constraint, this bound decreasing with increasing Aq. (We note 
that this lower bound is sensitive to the NLO corrections discussed in Sec. 2 
above.) Since larger Aq allows for larger ttiq in the coannihilation region, the 
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scattering cross section is a decreasing function of A . This is shown in Fig. 6 
where c^o_ p is plotted as a function of m 1( /2 for tan/3 = 40 and A = 2mi/2, 
4m 1/2 . 
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Figure 6. as a function of 171^2 for mSUGRA, fi >0, tari/9 = 40 and Ao = 2m^/2 

(upper curve), 4771J/2 (lower curve). 

We consider next /1 < 0. As discussed in 111, for low and intermediate tan (3, 
an accidental cancellation can occur in the heavy and light Higgs amplitudes 
in the coannihilation region which can greatly reduce a^o_ p . We investigate 
here what happens at larger tan/3, and what is the domain over which this 
cancellation occurs. In Fig. 7 we have plotted o~^o_ p in the large region, 
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for tan/3 = 6 (short dash), lO(solid), 20(dot-dash), and 25(dashed). One sees 
that the cross section dips sharply for tan/3 = 10, reaching a minimum at 
TO i/2 =725 GeV, and then rises. Similarly, for tan/3 = 20, the minimum 
occurs at mx/i = 830 GeV while for tan/3 = 25 at = 950 GeV. As a 
consequence, o~^o_ p wm below the sensitivity of planned future detectors 
for >450 GeV in a restricted region of tan/3, i.e. 

o-o_ p < 1 x 10 -10 for 450 GeV < m 1/2 < 1 TeV; 5 ~ tan/3 ~ 30; n < 0. (8) 

At the minima, the cross sections can become quite small, e.g. 1 x 10 _13 pb, 
without major fine tuning of parameters, corresponding to almost total can- 
cellation. Further, the widths of the minima at fixed tan (3 are fairly broad. 
While in this domain proposed detectors would not be able to observe Milky 
Way wimps, mSUGRA would imply that the squarks and gluinos then would 
lie above 1 TeV, but at masses that would still be accesible to the LHC. Also 
mSUGRA implies that this phenomena can occur only in a restricted range 
of tan /3, and for /i <0, so there would still be a number of cross checks of the 
theory. 




Figure 7. (°^o„ p ) for mSUGRA for fi < 0, A =1500 GeV for tan/3 = 6 (short dash), 10 

(solid), 20 (dot-dash), 25 (dashed). Note that the tan j3 = 6 curve terminates at low roi/2 
due to the Higgs mass constraint, and the other curves terminate at low rai/2 due to the 
b — » s-y constraint. 



In Sec. 2 it was pointed out that there was considerable uncertainty in the 
properties of the nucleon, particularly in the value of the av/y. In the above 
curves we have used a^N =65 MeV. Fig. 8 gives a comparison for tan (3 = 10 
of this choice (solid curve) with the parameters of E£l (dashed curve) where 
&ttN =45 MeV is used. One sees that the position of the minimum at mi/2 = 
725 GeV is shifted to 600 GeV, with similar shifts occuring for the other values 
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of tan/3. Thus the cancellations ocurring in this region are quite sensitive to 
the properties of the proton. 




450 500 550 600 650 700 750 800 

(GeV) 



2 

Figure 8. (o" s o_ p ) for mSUGRA for fj, < tan/3 = 10, a^N =65 GeV (solid), and the 
parameters of 12 (dashed) where ct^n =45 GeV was used. 



4 Nonuniversal SUGRA Models 

In most discussions of SUGRA models with nonuniversal soft breaking terms, 
the universality of the the soft breaking masses at Mq of the first two gen- 
erations of squarks and sleptons is maintained to suppress flavor changing 
neutral currents. However, one may allow both the Higgs masses and the 
third generation squark and slepton masses to become nonuniversal at Mg. 
One can parameterize this situation at Mq as follows: 

=m§(l + <yi); mj a = mjj(l + (9) 

m q L = m o(! + £3); m t l , =mo(l + Aj); m T 2 n , = mg(l + <5 5 ); 

TO fc« = m o(! + <*e); mjf = mo(l + £ 7 ). 

where qj, = (tt,&r,) squarks, II = (p T ,fi,) sleptons, etc. and too is the 
universal mass for the first two generations of squarks and sleptons. The 
6i are the deviations from universality (and if one were to impose SU(5) or 
SO(10) symmetry one would have <53=<54= 65, and Se=Sj.) In the following 
we limit the Si to obey: 

- 1 < Si < +1 (10) 

and maintain gauge coupling constant unification and gaugino mass unifica- 
tion at Mq, 
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While there are a large numbers of new parameters, one can get an un- 
derstanding of what effect they produce from the following. The neutralino 
Xi is a mixture of gaugino (mostly Bino) and higgsino parts: 

X° = aW 3 + (3B + 1 H 1 + SH 2 (11) 

The dominant spin independent c-oo-p cr °ss section is proportional to the 
interference between the gaugino and higgino amplitudes, and this interfer- 
ence is largely governed by the size of /j, 2 . As /Lt 2 decreases, the interference 
increases, and hence cr^o_ p increases. Radiative breaking of SU{2) x U(l) 
determines the value of fi 2 at the electroweak scale. To see the general nature 
of the effects of nonuniverality, we consider low and intermediate tan (3 where 
an analytic form exists for /i 2 (see e.g. Arnowitt and Nathu): 



t 2 - 1 



( 2 t 2 ' — 2 — ^ 3 + Si > ( 12 ) 



l + D St. 



universal parts + loop corrections. 



Here t = tan/3 and D = 1 - (m t /200GeVsin/3) 2 < 0.2 (Note that the Higgs 
and squark nonuniversalties enter coherrently, roughly in the combination 
S3 + 54 — S 2 ) We see from Eq.(12) that [i 2 is reduced, and hence er^?-p 
increased for S3, S4, Si < 0, 82 > 0, and /i 2 is increased for <5 3 , £4, Si > 0, 
S2 < 0. Thus one can get significantly larger cross sections in the nonuniversal 
models with the first choice of signs for the Si, and one can reduce the cross 
sections (though not by such a large amount) with the second choice. In 
general this would allow one to significantly lower the value of tan from the 
rquirement of tan (3 ~ 25 for mSUGRA to come within the range of current 
detectors. The matter is complicated, however, by the experimental Higgs 
mass constraint (that > 110) since theoretically the Higgs mass is small 
for low tan/3. We find, however, that one can reduce tan/3 to 7, and still 
maintain > 110 GeV. This is exhibited in Fig. 9, where the maximum 
cross section is plotted as a function of m^o . The nonuniversal cross sections 
can be a factor of 10 or more greater than the corresponding universal ones 
allowing for these much lower values of tan/3. 

As in mSUGRA, the minimum cross sections occur for the largest mi/2 
and smallest tan /3, and so they occur in the coannihilation region. We consider 
here the case where only the Higgs masses are nonuniversal i.e. Si t 2 7^ (the 
other S{ set to zero). Results then are similar to the mSUGRA case. For 
fi > we find 

a^o-p ~ 0.3(1) x 10" 10 pb; for m 1/2 < 1(0.6) TeV (13) 
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Figure 9. Maximum <r-o_ p for nonuniversal model for fi >0 with S3, 84, <5i < 0, 82 >0. 
Lower curve is for tan = 7 and the upper curve is for tan j3 = 12. 

For /i <0, one can again get a cancelation reducing the cross sections below 
1 x 10~ 10 pb for mi/2 > 450 GeV. As in mSUGRA these can produce sharp 
minima in the cross sections in the region tan/3 =10-25. 

5 D-Brane Models 

Recent advances in string theory has stimulated again the building of string 
inspired models. We consider here models based on Type II B oripntifolds 
where the full D = 10 space was compactified on a six torus T 6 u. These 
models can contain 9-branes and 5-branes which can be used to embed the 
Standard Model. We consider here a model in which SU(3)c x U(1)y is 
associated with one_set of 5-branes, 5i and SU{2)l is associated with a second 
intersecting set 52 □• Strings beginning and ending on 5i will have massless 
modes carrying the SU (3)c x U (l)y quantum numbers (i.e. the R quarks and 
R leptons), while strings starting on 52 and ending on 5i will have massless 
modes carrying the joint quantum numbers of the two branes (i.e. the quark, 
lepton and Higgs doublets). This then leads to the following soft breaking 
pattern at Me- 
nu = m 3 = -A = v / 3cos6'f ) 6ie" Ml m 3 /2 (14) 
m 2 = V3cosfl b (l-e?) 1/2 m 3/2 
where and rhi are the gaugino masses, and 

m\ 2 = (1 - 3/2 sin 2 9 b )m 2 3/2 for q L , l L , H x , H 2 (15) 
to 2 = (1 — 3 sin 2 8b)iriy 2 for ur, d,R, eR. 
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Thus the SU (2) doublets are all degenerate at Mq but are different from the 
singlets. We note Eq. (15) implies 9b <0.615 

This model was initially studied to examine its CP violating properties. 
However, it was subsequently seen that the experimental constraint of the 
electron electric dipole moment led to a serious fine tuning problem at Mq 
unless tan [3 ~ 3 — 5 El. Since we are interested here in larger tan /3, we will 
set the CP violating phases to zero. 

In general, Eq (15) shows that cr^o_ p is an increasing function of 9b since 
the squark and slepton masses decrease with increasing 9b . Thus the maximum 
cross sections will arise from large 9b and large tan/3. This is illustrated in 
Fig. 10, where o~^o_ p is plotted as a function of m^o for /i > for tan/3 =20, 
and 9b = 0.2. Thus we see that current detectors obeying the bound of Eq. 
(1) are sampling the parameter space for 



tan (3 ~ 20 



(16) 



We note that when tan /3 is close to its minimum value, m^o is also close to 
it's current LEP bound of m^o >37 GevElfi The minimum value of cr^o_ n 



3* 



Figure 10. <r~o_ 



mjjO (GeV) 



for D-brane model for fi > 0, 0^ = 0.2 and tan/3 = 20. The gap in the 



curve is due to excessive early universe annihilation through s-channel Z and h poles. 



will occur at low 9b, low tan/3, and large m 3 / 2 (i.e. large m^o). In the large 
m^o region, coannihilation can occur between the sleptons and the neutralino 
in a fashion similar to the SUGRA models, with the effective slepton mo 
parameter and effective neutralino m-i/a parameter being given by 

ttjq = (1 — 3 sin 2 9b) my 2 (17) 
mi/2 = V% cos 9b 0i m 3 / 2 
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Figure 11. Minimum <r-o_ for the D-brane model for a < and tan/3 = 6. 



Fig. 11 exhibits the minimum cross section for (i > as a function of the 
neutralino mass. One sees that 

cr^o_ p ~ 1 x 1(T 9 pb for n > (18) 

which is accessible to planned detectors. We note also that coannihilation is 
possible between the light chargino and neutralino. However, this occurs for 
only a very small region of parameter space. 

As in mSUGRA, a cancellation of matrix elements can occur for fi < 0, 
allowing for the cross sections to fall below the sensitivities of planned future 
detectors. This is exhibited in Fig. 12, where &^o_ p is plotted for tan/3 = 6 
(solid curve), 12 (dot-dash curve), and 20 (dashed curve). (The tan/3 = 6 
curve terminates at low m^o due to the constraint, while the higher tan /3 
curves terminate at low m^o due to the b — > sj constraint. The upper bound 
on m^o, corresponding to trig <1 TeV, arises from the Qh 2 constraint.) One 

sees that the cross section goes through a minimum at tan/3 = 12, though 
the expected rise at higher m^o does not appear since the parameter space 
terminates before this sets in. 

6 Conclusions 

We have examined here the neutralino proton cross section for a number 
of models possesing grand unification of the gauge coupling constants at 
Mq = 2 x 10 16 GeV. The mSUGRA model possesses universal supersym- 
metry soft breaking masses. In the nonuniversal models, this universality is 
maintained for the first two generations of squarks and sleptons and for the 
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Figure 12. o^o- for the D-brane model for fi < and tan j3 = 6 (solid), tan j3 = 12 
(dot-dash), and tan/3 = 15 (dashed). 

gaugino masses, but relaxed for the third generation and the Higgs scalar 
masses. In the D-brane models, universality is violated for both squark and 
slepton masses at Mq as well as for gaugino masses. Thus each model has a 
different pattern of soft breaking. Our analysis here also includes the recent 
experimental bounds on the Higgs mass of m h > 114 GeV 113 (which for GUT 
models holds for all tan/3), and the recent theoretical determination of the 
large tan /3 corrections to the NLO b — > S7 decay Lrtl, both of which produce 
significant effects on limiting the SUSY parameter space. Despite the phys- 
ical differences between the different models, the general picture resulting is 
somewhat similar. Thus current detectors obeying Eq(l) are sensitive to sig- 
nificant parts of the parameter space. For mSUGRA they are sampling the 
regions where tan/3 ~ 25. The nonuniversal models can have cross sections a 
factor of 10 or larger (with an appropriate choice of nonuniversalities) and so 
can sample the parameter space with tan/3 ~ 7. The D-brane models require 
tan/3 ~ 20. 

Coannihilation effects play a crucial role for large m^o in all the models, 
and for large tan/3, they are sensitive to the value of Aq. Large Aq leads to 
coannihilation corridors where too can get quite large, thus lowering the value 
of the <JyXi_ p cross section. For fj, > 0, the cross sections will generally still be 
accessible to planned future detectors obeying Eq. (2), i.e. 

a^o_ p ~ 1 x 10~ 10 pb for m 1/2 < 1 TeV, (i > (19) 

However, in all models, a special cancellation of the Higgs amplitudes can 
occur for fi < allowing the cross section to fall below the above bound when 
m i/2 ~ 450 GeV. For mSUGRA, these cancellations produce minima where 
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the cross section essentially vanishes for a range of m\/2 when 8 ~ tan f3 ~ 30, 
for mi/2 < 1 TeV (see Fig. 7) with similar results holding for the nonuniversal 

models. The cancellations for the D-brane models occur for 10 ~ tan/3 - 15. 
We note that at fixed tan (3 the cancellations can occur over a wide range of 
mi/j, ag. for mSUGRA tan/? = 10, er^o_ p < 10~ 10 pb for 400 GeV < m 1/2 < 
1 TeV Ej. In such regions of parameter space, dark matter detectors would 
not be able to observe Milky Way dark matter. However, these regions of 
parameter space would imply that gluinos and squarks lie above 1 TeV, but 
still should be accessible to the LHC if the parameter space is bounded by 
m i/2 < 1 TeV. Thus other experimental consequences of the models would be 
observable. 
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discussed in 

lower bound on cr^o_ p would be set by the spin dependent part of the 
cross section. Precisely when this would occur depends on the nuclei used 
in the target detector. 
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3, when the above cancelation is almost complete, the true 
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